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PHYSICAL  REVIEW  B,  in  press 

Quantum  Statistical  Theory  of  Laser-Stimulated  Surface 
Processes.  II.  Nonequilibrium  Rate  of  Laser  Excitation 
and  Multiphonon  Relaxation 

Jui-teng  Lin 
Laser  Physics  Branch 
Optical  Sciences  Division 
Naval  Research  Laboratory 
Washington,  D.C.  20375 

and 

Thomas  F.  George 
Department  of  Chemistry 
University  of  Rochester 
Rochester,  New  York  14627 

By  means  of  quantum  statistical  theory,  the  laser  excitation 
rate  and  multiphonon  relaxation  rate  are  calculated  for  species  ad¬ 
sorbed  on  cold  and  heated  surfaces.  The  complex  frequency  of  Bose 
operators  is  derived  from  a  microscopic  Hamiltonian  with  anharmonic 
mode-mode  coupling,  which  provides  the  origin  of  the  non-Arrhenius 
form  for  the  multiphonon  relaxation  rate.  The  energy  transfer  rate 
due  to  the  interference  effects  between  the  coherent  laser  excita¬ 
tion  and  multiphonon  coupling  is  calculated.  The  nonequilibrium 
energy  transfer  rate  due  to  the  energy  feedback  of  the  heated  sur¬ 
face  is  investigated  via  a  self-consistent  heat  diffusion  equation. 
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1 .  Introduction  | 

Laser-stimulated  surface  processes  (LSSP) ,  motivated  by  their 
potential  importance  in  industrial  applications  and  basic  research 
in  many  areas,  e.g.,  surface  chemistry,  materials  science  and  micro¬ 
electronic  processing,  are  receiving  increasing  attention.  LSSP 
may  be  grouped  into  the  selective  type,  e.g.,  excitation  of  a  spe- 
cific  active  mode,  adspecies  or  migrational  site,  and  the  non- 
selective  type  in  which  the  laser  radiation  energy  is  transformed 
into  heat. ^  These  selective  and  nonselective  aspects  of  LSSP  are 
characterized  not  only  by  the  pumping  rate,  as  compared  to  the 
multiphonon  relaxation  rate,  but  also  by  the  synergistic  effects 
of  the  coherent  laser  pumping  and  the  thermal  phonon- induced 
excitation.4 

Energy  transfer  processes,  resulting  from  laser  excitation  of 
an  adspecies  with  consequent  ntultiphonon  relaxation,  may  be  described 
by  the  rate  equation  for  the  active-mode  excitation  nA, 5 

nA<t)  *  WL  “  wRtnA(t)-nB]»  (i) 

where  W_  and  W_  are  the  laser  excitation  and  phonon-induced  relaxation  j 

"  A 

rate,  respectively,  of  the  active-mode  excitation,  whose  equilibrium 
value  nB  i®  given  by  a  Bose-Einstein  distribution.  The  above  rate 
equation  can  be  derived  from  a  microscopic  Hamiltonian,  where  the 
Markovian  approximation  is  assumed  such  that  the  "feedback"  energy 
from  the  heated  surface,  characterized  by  the  temperature-dependent 
bath-mode  occupation  nB(T(t)],  is  completely  eliminated,  i.e.,  an 
equilibrium  statistics  is  used  to  treat  the  laser-induced  transient 
phenomena.  Furthermore,  a  constant  multiphonon  relaxation  rate. 
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WR,  is  assumed  in  the  rate  equation. 

Energy  transfer  processes  may  also  be  pictured  by  a  master 
2  4 

equation/  ' 

K  -  I'vw.'  •  «J> 

in 

where  PR  is  the  energy  population  of  state  n  ,  and  Wmw  is  a  total 
transition  rate  for  a  transition  from  m  to  n  .  This  rate/ 
induced  by  both  the  coherent  laser  field  and  the  multiphonon  relaxa¬ 
tion,  cannot  in  general  be  separated  into  the  sum  (total)  ■ 

Wmn (laser)  +  W^n (phonon) ,  particularly  for  the  situation  of  a 
strong  phonon  coupling  and/or  high-power  laser  excitation  processes. 
Furthermore/  the  excitation  of  the  phonon  modes  can  result  in  a 
hot  surface  in  which  the  feedback  energy  of  the  thermal  phonons 
significantly  affects  the  overall  transition  rate. 

Here,  we  shall  develop  a  quantum  statistical  theory  for  the 
total  transition  rate  which  includes  the  interference  effects 
between  the  coherent  laser  field  and  the  incoherent  multiphonon 
coupling.  In  Section  2,  a  microscopic  Hamiltonian,  including  the 
enharmonic  mode-mode  coupling  which  provides  the  origin  of  a  non- 
Arrhenius  transition  rate  is  investigated.  The  equilibrium  transi¬ 
tion  rate  for  a  cold  surface  and  the  nonequilibrium  transition  rate 
for  a  heated  surface  with  energy  feedback  effects  are  studied  in 
Sections  3  and  4 ,  respectively.  Concluding  remarks  are  given  in 
Section  5. 


2.  Microscopic  Hamiltonian  with  Anharmonlc  Couplirv 


We  consider  a  system  consisting  of  species  adsorbed  on  a 

solid  surface  and  subject  to  infrared  radiation,  with  low  intensity 

for  cold-surface  processes  or  high  intensity  for  heated-surface 

processes.  The  spectrum  of  the  adspecies/surface  system  may  be 

partitioned  into  a  high-frequency  mode  A,  which  is  the  active 

vibrational  mode  of  the  adspecies,  and  a  low-frequency  region 

comprised  of  the  bath  modes  B  and  C.  Due  to  the  energy  gaps  among 

the  A,  B  and  C  modes,  the  photon  energy  is  selectively  transferred 

to  the  A  mode  with  subsequent  relaxation  to  the  B  and  C  modes  via 

mainly  multiphonon  processes.  The  total  vibrational  Hamiltonian 

2 

describing  the  above  mentioned  system  may  be  expressed  as 

H  *  SA  ♦  “b  *  HC  *  HAC  +  he  ♦  H»S  -  '  131 

where  Hi  (i»A,B,C)  are  the  unperturbed  Hamiltonians  of  the  A,  B 
and  C  modes,  respectively,  while  H^.,  and  H^  are  the  interac¬ 
tion  Hamiltonians  among  these  modes  [we  note  that  H^  can  include 
the  enharmonic  coupling  among  the  bath  modes  (B+C) ] .  H^,  is  the 
active  mode  and  laser  field  interaction  Hamiltonian  in  which  a 
classical  field  is  assumed,  a  schematic  diagram  of  the  laser  exci¬ 
tation  and  multiphonon  relaxation  processes  is  shewn  In  Figure  1(A)  along  with 
a  sketch  of  the  frequency  spectrun  of  the  system  in  Figure  1(B) .  The  eval¬ 
uation  of  the  pumping  and  relaxation  rates  will  be  based  on  the 

above  total  Hamiltonian. 

To  treat  the  many-body  problem  due  to  the  multimode  couplings 
one  may  use  a  Markovian  approximation  in  which  a  fast  characteristic 
time  of  the  bath  modes  is  assumed,  or  use  a  Wigner-Weisskopf  approxi 
mation  in  which  a  single  pole  is  assumed  in  the  Laplace  transform 
of  the  relevant  operators.6  With  either  of  these  techniques,  the 
effects  of  the  multiphonon  coupling  terms  H^c  and  HfiC  may  be  ef¬ 
fectively  absorbed  into  the  " C-mode- perturbed"  Bose  operators 


a  and  b  for  tha  A  and  B  modes,  respectively,  resulting  in  complex 
frequencies  whose  imaginary  parts  describe  the  energy  relaxation 
of  the  A  and  B  modes  into  the  C  modes.  Mathematically  this  can 
be  expressed  in  terms  of  the  Heisenberg  equations  of  motion  for 
a  and  b, 

i-ni  -  1SQ0  +  [a.H^+H^],  (4. a) 

ilibj  ■  Hi  Q  j+  Ibj,HAB], 

where,  within  the  Markovian  approximation,  the  effective  ( complex) 

frequencies  ftg  and  are  eigenfrequencies  of  the  C-mode-perturbed 

2  7 

Hamiltonians  HA  +  HAC  and  Hg  +  H^,  respectively,  given  by  ' 

fl0  “  “o  +  6u0  “  2e*nA  +  iY0/2,  (5. a) 

**  Wj  f  iYj/2<  (S.b) 

Here  u>Q  and  are  the  eigenfrequencies  of  HA  (with  anharmonicity 
e*  ud  excitation  nA)  and  H0  (j-th  mode),  respectively,  and  yQ(&wQ) 
and  y j  (6u>j )  are  the  C-mode-induced  "level  broadening"  (frequency 
shift) ,  including  both  (energy)  and  T2  (phase)  relaxation,  of 
the  A  and  B  modes,  respectively.  ftQ  and  £)j  may  be  simply  phenomen¬ 
ologically  chosen  or  rigorously  derived  from  a  microscopic  Hamiltonian. 
From  the  latter  method  the  damping  factor,  y^  is  found  to  be  related 
to  the  density  Of  states  of  the  C  modes  and  the  mode-mode  coupling 

strength,  K,  by  y0  "  2ir|K|  p,  for  single-phonon  coupling,  and  a 

2 

convolution  form  for  multiphonon  coupling. 

We  now  rewrite  the  total  Hamiltonian  [Eq.(3)]  as 


-6- 


H  -  HI  +  V, 


(6. a) 


HA  »  {H4+H- _)  +  (HB+H_  ) , 


A  AC 


B  BC 


(6.b) 


v“hab  +  HAF' 


(6.c) 


in  which  the  Bose  operators  in  an  interaction  picture  are  given 
in  terms  of  the  complex  frequencies ,  i.e.,  aI(t)  -  exp(-iH^tA)* 
a-expUH^t^)  has  the  expression 


and  similarly. 


aT(t)  •  a(0)  exp(-i  a.  t) , 


^(t)  »  bj  (O)exp(-iO*  t) . 


(7. a) 


( 7 .  b) 


We  note  that  these  Bose  operators,  with  complex  frequencies  in¬ 
duced  by  the  coupling  to  the  C  modes,  provide  us  with  a  novel 
feature  of  the  transition  rate  which  turns  out  to  be  non-Arrhenius. 


3.  Rates  of  Laser  Excitation  and  Multiphonon  Relaxation  for  a  Cold 
Surface 

We  shall  now  develop  a  quantum  statistical  transition  rate 
based  on  the  Hamiltonian  given  by  Eq. (6)  and  the  corresponding 
Bose  operators  in  the  interaction  picture.  Because  of  the  time- 
dependent  laser  interaction  Hamiltonian,  H  ( t) ,  the  transition 

Ar 

rates  are,  in  general,  time-dependent.  In  particular,  for  the 
situation  of  high-power  laser  excitation  with  strong  phonon 
coupling,  the  temperature  of  the  substrate  and  thus  the  occupation 
numbers  increase  during  the  pumping  process  of  the  active  mode. 

The  nonequilibrium  transition  rate  will  be  discussed  in  the  next 
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section.  We  shall  first  study  the  equilibrium  processes  for  a 
cold  surface,  where  a  constant  Bose-Einstein  distribution  is 
used  for  calculating  the  ensemble-averaged  quantities. 

Employing  Fermi's  Golden  Rule,  the  total  energy  transfer 
rate  of  the  active  mode  pumped  by  coherent  radiation  with  con- 

Q 

sequent  multiphonon  relaxation  governed  may  be  written  as 

,  •  I 

W  -  (~)R e/  dt «V(t)V(0)»,  (8) 

1i 2  0 

where  <<•••>>  denotes  an  ensemble  average  over  both  the  active 
mode  and  bath  mode  coordinates.  This  rate  is  characterized  by 
the  correlation  function  of  V(t)  in  the  interaction  picture,  i.e., 
V^t)  ■  exp(-iH^)V(t)exp(iH^tA) .  Knowing  that  the  active  mode  is 
perturbed  by  two  interactions,  the  coherent  radiation  and  the 
phonon-induced  relaxation,  we  may  decompose  the  total  transfer 
rate  into  three  components?  the  laser  excitation  rate  WL,  the 
multiphonon  relaxation  rata  WR  and  the  interference  term  W^: 


where 


WR  * 
MRL  * 


(t)  H 


AF 


(t)H 


AB 


(t)H 


AB 


(9) 

(0)»,  (10. a) 

(0) >>,  (10 .b) 

(0)  (io.c) 


We  first  calculate  the  laser  excitation  rate  W^.  For  a 
classical  radiation  field  E(t)  ■  EgCosw^t,  we  have 
HAp(t)  -  Vg  (a*+  a)cosu>Lt,  where  Vg  is  proportional  to  the  electric 
field  Eg  and  the  derivative  of  the  active  dipole  moment.  Combining 
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the  expression  in  E^.(7.a)  end  its  Hermitian  conjugate,  we  can 
write  HAJ1(t)  in  the  interaction  representation,  using  the  rotating- 
wave  approximation ,  as 

HAF(t)  *  (V0/2)  (a+(0)exp[i(fl0-u»L)t]+  a(0)exp[-i(Q*-«L)t]).  (11) 


Substituting  Eq.(ll)  into  Eq. (10. a)  and  working  out  the  integration, 
we  obtain 


W, 


2  ( 


(12) 


where  nA  *  [exp  (fta)g/kT)  -1]  is  the  initial  occupation  number  of  the 
active  mode  and  A  ■  u0~ul  is  the  detuning.  In  obtaining  Eq.(12), 
a  low  excitation  was  assumed  such  that  the  time-dependent  detuning 
A(t)  -  A  -2e*nA(t)  is  approximated  by  A.  For  high  excitation  or 
large  enharmonic ity,  a  time-dependent  excitation  of  the  active  mode, 
nA(t) ,  must  be  solved  in  order  to  obtain  the  laser  escitation  rate. 
An  exact  solution  in  this  case  requires  a  numerical  integration 
which  will  not  be  discussed  here.  However,  using  an  iterative  pro¬ 
cedure,  we  find  that  the  first-order  correction  to  WL  is  given  by 
the  same  Lorentzian  form  as  that  of  Eq.(12)  with  a  modified  detuning 


A(1)  *  A  -  2e*(v0/2R)2/[A2+(Y0/2)2]. 


To  calculate  the  multiphonon  relaxation  rate  from  Eq.(lO.b), 
the  interaction  Hamiltonian  is  required.  We  consider  a  multi- 

Q 

phonon  coupling  Hamiltonian  of  the  form7 

HAB  "  9,(a++a)  [exp (-lag)-!.] ,  (la) 

where  g  is  a  coupling  strength  and  Q  is  the  normal  coordinate  of 


the  surface  atom  coupled  to  the  adspecies.  Rewriting  ^  in  the 

interaction  representation,  substituting  it  into  Eq.(lO.b)  and 

using  the  cumulant  expansion  for  the  correlation  function, 

2 

<<exp(i£)>>  sexp<<-£  /2>>,  we  obtain 

2  •  r  -in*t  iarttT 

**R  -  2  If  I  *•/  dt[(5A+l).  *nAe  J- 


r  G.(t)  +  G_(t) 

(e  '(«  +  )-(eTt.l)e  +1 

]. 

(14) 

•  2  _  ifl(w.jt 

G+(t)  ■  /  dUjp(<4ij)8^(<Uj)n(Uj)e  J  , 

(14. a) 

<*>  _  -ift*(w.)t 

G_(t)  -  /  du)jP(ia^)  [H(Wj)+l]e  •  3  , 

(14. b) 

G  -  G+(0)  +  G_(0)  , 

(14. c) 

where  the  complex  frequencies  and  are  given  in  Eq.(5)r 

6  (ojj)  is  a  quantization  factor  resulting  from  Q  »  E[6(ujj)/o]  (b^+bj) 
and  p(uj)  is  the  density  of  states  of  the  B  modes  with  a  mean 
damping  factor  YjxT  and  occupation  number  n(aij)  for  a  continuous 
spectrum. 

The  above  integral  expression  describing  the  multiphonon 
relaxation  rate  of  the  active  mode  for  an  arbitrary  value  of  G  can¬ 
not  be  solved  analytically.  However,  for  small  G,  we  may  expand 
the  exponent  of  the  integrand,^0 


explG, (t)+G_(t) ]  -  [G  (t)+G_(t) ]m, 
+  ra-0 


(15) 


to  obtain  the  single-phonon  (m*l)  and  multiphonon  relaxation  rates. 


For  an  Einstein  spectrum,  p(u>4)  «  5  (u)j-u>Ei ,  we  obtain  the  result 
_  2  -  _  N 

WR  -  2|J|  e.<*(2nA+l)  Z  L(Dm)I(x),  (16. a) 

m*0 


L<Dm)  -  V<D*«i>.  (16  .b) 

wher.  I(x) ,  with  x  ■  26 tn(n+l) ]2^2,  and  6  =  0(uE) ,  is  a  modified 

Bessel  function,  Dffl  is  the  m-phonon  detuning  given  by  Dm  *  a»0-ma»E 

and  the  overall  damping  factor  is  r+  -  [(2-G)7+Ygl/2.  In  writing 

Eq. (16) ,  we  kept  only  the  dominant  term  and  discarded  the  off- 

resonance  terms  involving  the  Lorentzian  functions  L(u»q)  and 

L«J>,  -  u>Q-Hnu>E.  We  note  that  the  multiphonon  relaxation  rate 

for  the  small  G  case  is  Arrhenius -like  due  to  the  exponent  exp(-G) 

for  a  high- temperature  limit  in  which  G«T.  However,  WR  is  a 

2 

Lorentzian  in  the  low- temperature  limit  where  G  s  B  u»E  becomes 
independent  of  temperature. 

We  now  investigate  the  large  G  value  case.  Instead  of  taking 
the  expansion  of  exp(-G)  as  we  did  for  small  G,  we  employ  the 
steepest  descent  method  in  which  the  short-time  value  of  G  is  the 
dominant  contribution.11  The  result  for  large  G,  in  an  Einstein 
spectrum,  is  obtained  from  Eq.(14)  as 

WR  -  2||12(EA+1)L(DS) ,  (17. a) 

L(DS)  -  ( Yg/2) / (Dg+ (Yg/2) ^] ,  (17 ,b) 

2 

where  the  detuning  is  defined  by  Dg  »  u>0  -  8  wE,  and  again  we  have 

discarded  the  off-resonance  terms  involving  Mu>q)  and  L(Dg), 

4*  2 

Dg  ■  +  6  u>E<  We  note  that  the  above  non-Arrhenius  form  for  the 

multiphonon  relaxation  rate  is  significantly  different  from  that 
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12 

of  previous  work,  and  the  origin  of  this  I-*' cant zian- type  relaxa¬ 
tion  rate  is  the  enharmonic  couplings  of  the  A  and  B  inodes  to  the 
C  modes  which  cause  the  damping  of  the  (A+B)  modes.  In  the  limit 
of  very  weak  C-mode  coupling,  i.e.,  r  ♦  0,  the  Debye-Waller  factor 
(proportional  to  G)  contained  in  the  exponent  exp(-G)  would  not 

be  completely  cancelled  out  as  in  the  r+yo  case,  and  Eq.  (14) 

12 

results  in  a  usual  Arrhenius  form. 

Finally,  the  interference  term  given  by  Eq.(lO.c)  may  be 
calculated  by  using  the  interaction  representation  for  H^,  [eq.(ll)] 
and  [eq.(13)].  The  result  is 


WRL  "  (2nA+l)L(A)  (l-ft'G/2) , 

■ft 

L(A)  -  (Y0/2)/[A2+(Y0/2)2J, 


(18. a) 
(18. b) 


for  arbitrary  values  of  G,  where  A  »  wQ  -  <*»L  is  a  laser  detuning. 
Note  that  the  above  interference  term  provides  an  extra  contribu¬ 
tion  to  the  excitation  of  the  active  mode,  which  is  "constructive" 
for  l>exp(-G/2),  and  "destructive"  for  l<exp(-G/2).  This  con¬ 
tribution  will  be  significant  (i.e.,  the  total  energy  transfer 
rate  cannot  be  simply  decomposed  into  two  parts  w  ■  +  WR) 

particularly  when  the  factor  (VQg)L(d)  in  W^  is  comparable  to 
the  factor  |g|2L(D)  in  WR. 


4 .  Nonequilibrium  Energy  Transfer  Rate  for  a  Heated  Surface 

In  the  previous  section,  the  energy  transfer  rates  were 
calculated  for  a  system  with  a  cold  surface,  such  that  the  Bose 
functions  nA  and  n(w)  maintained  the  equilibrium  values  without 
significant  increase  during  a  low-power  laser  excitation.  However, 


-12- 


£or  a  system  initially  at  room  temperature  (TQ)  exposed  to  high- 
power  laser  radiation,  a  change  of  temperature  from  TQ  « fia>E  to 
T,  HuE<<kT  <‘Huq »  is  significant,  and  the  constant  Bose  functions 
must  be  replaced  by  their  transient  values,  namely  nA(t)  and 
n(u>/ 1) .  The  nonequilibrium  energy  transfer  rate  characterized  by 
the  transient  Bose  functions,  which  in  turn  are  given  by  the  instant 
temperature  of  the  system,  may  be  obtained  by  solving  a  self-consistent 
heat  diffusion  equation3'13 


3T 

St 


7- (D7T)  +  WR(T, t) , 


(19) 


where  T  is  the  laser-induced  transient  temperature  of  the  system 
with  diffusivity  D  and  heat  capacity  Cv,  and  heat  (energy)  diffuses 
according  to  the  gradient,  7,  defined  in  the  direction  of  energy 
flow.  0  is  a  coverage  factor  given  corresponding  to  the  number  of 
active  sites  within  the  interaction  region.  The  above  equation  is 
highly  nonlinear  because  of  the  complicated  temperature  dependence 
of  the  multiphonon  relaxation  rate  HR[nA(T)  ,n(oi,T)  ] ,  and,  in  general, 
the  diffusivity  is  also  temperature  dependent,  e.g.,  D  »  DQ  +  D  T, 
for  most  metal  substrates. 

An  exact  solution  to  Eq . ( 19 )  requires  a  numerical  integration 

which  is  not  presented  here.  Instead,  to  investigate  the  effects 

of  the  nonequilibrium  Bose  function  on  the  energy  transfer,  we 

study  the  tractable  case  with  T  »  Tq  +  At,  where  TQ  is  an  initial 

0  0 

temperature  and  A  -  0WR/Cv,  with  WR  being  the  initial  relaxation 
rate  given  by  Eq.(17.a).  Using  this  linearly  time-dependent  tem¬ 
perature,  we  are  able  to  calculate  the  energy  transfer  rates  including 
the  energy  feedback  effects  from  the  laser-heated  surface.  For 


a  large  G  value  with  an  intermediate  temperature  RwL<<kT<t&«>0, 
we  obtain  the  multiphonon  relaxation  rate  in  the  presence  of  sur¬ 
face  heating  up  to  the  first-order  as 


-•( 


o  e  r^-D2 

p  0  V  m 


1  +  (c  t'  * 


r«>. 


!>4 


(21) 


where  B  *  £  -  t  (2— G)  y+YqI/2  and  Dm  *  “0  “  ®«E  i*  the 

detuning.  The  laser  excitation  rate  with  energy  feedback  from  the 
thermal  phonons  is  found  to  be 


W 

i 

where  6y 
effects. 


(1) 


1 1  V°  I  2  f  n 
zi-s-i 


vo . 2  f-  ,V  W 
'75 


A '  A2  + (  y  - 6 Y )  2/ 4 
0 


)  + 


<V2) 


-5 - J 

A2+(yq/2)2 


* 

r  * 


(22) 


20W°/Cv  is  a  correction  factor  caused  by  the  heating 


Finally,  the  interference  term  is  given  as 

WRL)  “  WRL[1  +  6y  (23) 

Comparing  the  above  transfer  rates  for  a  heated  surface 
[Eqs. (21)-(23) ]  with  those  of  a  cold  surface  [Eqs. (12) , (17)  and 
(18)],  we  see  that  as  the  surface  temperature  increases,  both  the 
laser  excitation  rate  W_  and  the  multiphonon  relaxation  rates 

ii  • 

WR  and  are  enhanced.  Therefore,  there  is  a  competition  between 
the  feedback-enhanced  multiphonon  relaxation  and  the  laser  excitation 
rate  when  the  surface  is  significantly  heated.  For  large  photon 
energy  stored  in  the  active  mode,  which  is  the  essential  component 
of  selective  type  processes,  a  cold  surface  is  required  where  the 
coherent  excitation  will  not  be  smeared  out  by  feedback  (incoherent) 
effects.  However,  for  a  nonselective-type  excitation,  the  feedback 
energy  provides  an  extra  contribution,  in  addition  to  the  laser 
source,  to  either  or  W ^ . 


5.  Conclusion 

From  a  microscopic  C-raode-perturbed  Hamiltonian,  the  complex 
frequencies  of  the  Bose  operators  are  derived  and  shown  to  be  the 
origin  of  the  non-Arrhenius  form  for  the  multiphonon  relaxation 
rate  (Eq.(5)].  There  is  an  interference  term  accompanying  the 
coherent  laser  excitation  and  the  multiphonon  relaxation,  and  the 
contribution  of  this  term  may  be  significant  for  strong  phonon 
coupling  and/or  strong  laser  excitation  (Eq.(18)].  For  high- 
power  laser  excitation,  there  is  a  nonequilibrium  energy  transfer 
rate  which  is  characterized  by  the  transient  temperature  of  the 
system.  The  energy  feedback  from  a  heated  surface  may  be  described 
by  a  self-consistent  heat  diffusion  equation  [Eq.(19)].  For  a 
simple  example  with  temperature  linearly  increasing  in  time,  we 
have  shown  the  enhancement  of  both  the  laser  excitation  rate  and 
the  multiphonon  relaxation  rate  [Eqs. (21) -(23) ] . 

f 
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Figure  1.  (A)  Schematic  diagram  for  the  laser-excited  adspecies/surface 
system  in  terms  of  the  laser  pumping  rata  (WL) ,  the 
multiphonon  relaxation  rate  (WR)  and  the  damping  factors 
Yq  and  y j  for  the  active  and  bath  modes,  respectively. 

(B)  Sketch  of  the  density  of  states  of  the  system  as  a 

function  of  frequency.  By  the  concept  of  the  "energy-gap 
law”,  the  A-C  couplings  are  weaker  than  the  A-B  couplings. 
Furthermore,  since  there  is  no  overlap  between  the  A  and 
B  modes,  a  single-phonon  process  is  forbidden  so  that 
the  relaxation  rate  depends  only  on  multiphonon  processes. 
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